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Abstract 

The meson-baryon coupled channel unitary approach with the local hidden gauge formalism is 
extended to the hidden beauty sector. A few narrow N* and A* resonances around 11 GeV are pre- 
dicted as dynamically generated states from the interactions of heavy beauty mesons and baryons. 
Production cross sections of these predicted resonances in pp and ep collisions are estimated as a 



guide for the possible experimental search at relevant facilities. 
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I. INTRODUCTION 



In conventional quark models, all established baryons are ascribed into simple 3-quark 
(qqq) configurations The excited baryon states are described as excitation of individual 
constituent quarks, similar to the cases for atomic and nuclear excitations. However, unlike 
atomic and nuclear excitations, the typical hadronic excitation energies are comparable 
with constituent quark masses. Hence to drag out a qq pair from gluon field could be a 
new excitation mechanism besides the conventional orbital excitation of original constituent 
quarks. Some baryon resonances are proposed to be meson-baryon dynamically generated 
states or states with large (qqqqq) components }9l4ll|. A difficulty to pin down the 

nature of these baryon resonances is that the predicted states from various models are around 
the same energy region and there are always some adjustable ingredients in each model to 
fit the experimental data. A typical example is iV*(1535) which has large couplings to 
the strangeness. In the 3-quark (qqq) configurations, it is described as the orbital angular 
momentum L = 1 excitation of a quark. But phenomenological studies suggest that it may 



be a quasi-bound state of KT, system [l2Nl4j. or a hidden strangeness 5-quark state |lOl. Il5|. 
In order to clearly demonstrate the new excitation mechanism with some of its corresponding 
states, in Ref . 16j| , the meson-baryon coupled channel unitary approach with the local hidden 



gauge formalism was performed for the hidden charm sector and several narrow iV* and 
A* resonances with hidden charm were predicted to exist. If found experimentally, these 
resonances definitely could not be described as three constituent quark states. Here, we 
extend the study to the hidden beauty sector. Some super-heavy iV* and A* resonances 
with hidden beauty are predicted to exist, with mass around 11 GeV and width smaller 
than 10 MeV. If these resonances would be experimentally confirmed, they should be part 
of the heaviest super-heavy island of N* and A* state. As a guild to the future experimental 
search for these new predicted states, their production cross sections in pp and ep collisions 
are estimated. 

In the next section, we present the formalism and ingredients for the study of interactions 
between heavy beauty meson and baryon with the Valencia approach, and give some detailed 
discussion on the intermediate meson-baryon loop G functions. In section 1111} our numerical 
results for the masses and widths of the predicted super-heavy N* and A* states are given, 
followed by a discussion. In section HVl effects of momentum dependent terms in the effective 
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potential are investigated. In section [V] the calculation about production of these predicted 
states from pp and ep collisions is presented. Finally, a short summary is given in the last 
section. 



II. FORMALISM FOR MESON-BARYON INTERACTION 



We follow the recent work of Ref. 



16] on the interactions between charmed mesons 



and baryons, and replace charm quark by beauty quark. The PB — > PB and VB — > 
VB interactions by exchanging a vector meson are considered, as shown by the Feynman 
diagrams in Fig. [TJ 




b 2 b; 




(a) 



(b) 



FIG. 1: Feynman diagrams for the pseudoscalar-baryon (a) or vector-baryon (b) interaction via 
the exchange of a vector meson (Pi, P2 are B°, B + or B®, and V\, V2 are B *, B + * or B®*, and 
Bi, B2 are A&, H&, E' b or fib, and V* is p, K* , <f> or uS). 



The effective Lagrangians for the interactions involved are 171 ] : 



Cyyy = ig(V»[V, d^]) 

Cppv = -ig{V»[P,d„P\) 

Cbbv = gdB^V^B}) + {B^B){V»)) 



(1) 



where P and V stand for pseudoscalar and vector mesons of the 16-plet of SU(4), respectively. 



Using the same approach of Ref. 



16], we keep only the 7 component of Eq.flT]), while 



the three momentum versus the mass of the meson can be neglected under the low energy 
approximation. Similarly, the q 2 /M v term in the vector meson propagator is neglected so 
that the propagator is approximately equal to g^ u /My. Note when we consider transitions 



from heavy mesons to light ones later on, we perform the exact calculation without such ap- 
proximation. Then with g = My/2f the transition potentials corresponding to the diagrams 
of Fig. [TJ are given by 



ab(P 1 B 1 ^P 2 B 2 ) 



^(E Pl + Ep 2 ), 



Vn 



ab(ViB 1 ^V 2 B 2 ) 



ah (T? 



Ev 2 )ti ■ e*2, 



(2) 
(3) 



where the a, b stand for different channels of Pi(Vi)Bi and ^2(^2)52, respectively. The E is 
the energy of corresponding particle. The e is the polarization vector of the initial or final 
vector. And the e\ 2 component is neglected consistently with taking p/My ~ 0, with p the 
momentum of the vector meson. Here we only change the charm quark to beauty quark, 
so the C a b coefficients are exactly the same as those in Ref. so that there are only two 
cases, (I, S) = (1/2, 0) and (0, -1), which have attractive potentials. We list the values of 
the C a b coefficients for PB — > PB for these two cases in Table I and Table II, respectively. 



TABLE I: Coefficients C ab in Eq. © for (I, S) = (1/2,0) 
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TABLE II: Coefficients C ab in Eq. © for (J, S) = (0, -1) 
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With the transition potential, the coupled-channel scattering matrix can be obtained by 



solving 
of Refs. 



;he coupled-channel Bethe-Salpeter equation in the on-shell factorization approach 



5] 



T = [1 — VG]- l V 
4 



(4) 



with G being the loop function of a meson (P), or a vector (V), and a baryon (B). The el • e*2 
factor of Eq. ([3]) factorizes out also in T. 

For the G loop function, there are usually two ways to regularize it. The first one is using 
dimensional regularization by means of the formula 

d A q 1 1 



G = i2M B 



(2tt) 4 (P-q) 2 -M 2 B +ie q 2 -M 2 p + ie'' 



2M B ( Ml M 2 -M 2 +s Mg 

- i 67r 2^ + m ^ + 2s M% 

+4= [ln(s - (M 2 B - M 2 P ) + 2q^7s) + \n(s + (M 2 B - M 2 P ) + 2q<Js) 

-\n(-s - (Ml - M 2 P ) + 2q^~s) - ln(-s + (M 2 B - M 2 P ) + 2qy/s)] } , (5) 

where q is the four-momentum of the meson, P is the total four-momentum of the meson 
and the baryon, s = P 2 , q denotes the three momentum of the meson or baryon in the 
center of mass frame, /i is a regularization scale, which we put 1000 MeV here. Changes in 
the scale are reabsorbed in the subtraction constant a M to make results scale independent. 
a M is of the order of —2, which is the natural value of the subtraction constant [18j. When 
we look for poles in the second Riemann sheet, we should change q to — q when a/s is above 
the threshold in Eq.© [id]. 

The second way to regularize the G loop function is by putting a cutoff in the three- 
momentum: 

d 4 q 1 1 



G = %2M B 



(2tt) 4 (P - q) 2 - M 2 B + ieq 2 - M 2 p + is 
A q 2 dq 2M B (u P + u B ) 

4-7T 2 up u B (s — (ujp + uj b ) 2 + ie) 



(6) 



where up = \Jq 2 + Mp, u B = -y/g 2 + M|, and A is the cutoff parameter in the three- 
momentum of the function loop. 

Here we give some detailed discussion on these two types of G function. Firstly the free 
parameters are a M in Eq.© and A in Eq.©. The value of A is around 0.8 GeV, which are 
within the natural range for effective theories [5j. Then the a M parameter is determined by 
requiring that the two G functions from Eq.© and Eq.©take the same value at threshold. 
This value also leads to similar shape near threshold for the two G functions as shown 
in Figj2j In FigJ21 the real part and imaginary part of two G functions vs the energy 
difference between the center mass energy and the corresponding threshold for i^S, DH C 



and BUb channels are demonstrated. In Table UTTj the parameters for different G functions 
and channels are listed. While the imaginary parts of two G functions are exactly the same, 
there are some differences for the real parts of two G functions and the differences become 
bigger for heavier channels. For the same A value, the magnitude of a M depends on the 
threshold of channels and gets bigger for heavier channels. One point worth mentioning 
is that for the BUb channel the real part of the G function given by Eq.([3]) is larger than 
zero for energies more than 50 MeV below the threshold as shown in FigfSJ As we know, 
if the interaction is repulsive potential, i.e., the value of the potential V is positive, there 
should be no bound state. However, when the real part of G function is also positive below 
the threshold, the pole can still be found in the model T matrix with a repulsive potential. 
These poles far below threshold are beyond the valid region of the model approximation 
and should be discarded. Since varying the G function in a reasonable range does not 
influence our conclusion qualitatively, we present our numerical results in the dimensional 
regularization scheme with a M = —3.71, in this paper. 



-100 -75 -50 -25 25 50 75 100 -100 -75 -50 -25 25 50 75 100 




Mem. - Mthreshold(MeV) Mem. - Mthreshold(MeV) 



FIG. 2: The real part (left) and imaginary part (right) of two G functions vs the energy difference 
between the CM. energy and the threshold energy. The solid lines are for Eq.©, and dashed lines 
are for Eq.([S|). The thickest lines are for BY^b channel, the thinnest ones are for KT, channel, and 
middle ones are for DTi c channel. The used parameters are listed in the Table iriTl with A = 0.8GeV^. 

With the potential and G function fixed, the unitary T amplitude can be obtained by 
Eq.flU). The poles in the T matrix are looked for in the complex plane of y/s. Those appearing 
in the first Riemann sheet below threshold are considered as bound states whereas those 
located in the second Riemann sheet and above the threshold of some channel are identified 
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TABLE III: The parameters for two types of G functions in the cases of KT,, Z)S C and BT,^ 
interactions, with for Eq.(|5]) and A for Eq.([6]). The listed a M and A(GeV) give the same value 
of two G functions at the corresponding threshold. 





Threshold(GeV) 






« 7 , 






A(GeV) 




0.7 


0.8 


0.9 


1.0 


1.1 


BZ b 


11.087 


-3.679 


-3.715 


-3.751 


-3.786 


-3.822 


DT, C 


4.231 


-2.196 


-2.283 


-2.369 


-2.453 


-2.536 




1.688 


-1.297 


-1.463 


-1.619 


-1.766 


-1.905 



as resonances. As previously discussed, the poles will be kept only when the real part of 
Eq.flHJ) is negative. 

From the T matrix for the PB — > PB and VB — » VB coupled-channel systems, we 
can find the pole positions zr. Six poles are found in the real axes below corresponding 
thresholds and therefore they are bound states. For these cases the coupling constants are 
obtained from the amplitudes in the real axis. These amplitudes behave close to the pole 
as: 



9a9b 



(7) 



y/~S- Z R 

We can use the residue of T aa to determine the value of g a , except for a global phase. Then, 
the other couplings are derived from 

/ 9aT a b \ 



9b 



lim 



^j~s^>z R T aa 

III. NUMERICAL RESULTS FOR THE SUPER-HEAVY N* AND A* 



(8) 



Firstly, we discuss the (I, S) = (1/2, 0) sector. There are two channels, SS^ and BAb. 
The masses of these particles are taken from [l], mg = 5.279 GeV, itlb* = 5.325 GeV, 
m-£ b = 5.808 GeV and m\ b = 5.620 GeV. With the approach outlined in the last section, the 
obtained pole positions zr and coupling constants g a are listed in Tables IIVI for PB — > PB 
and VB — > VB. Because these poles are bound states for each channel, they have zero 
width when neglecting transitions mediated by t-channel exchange of heavy beauty mesons. 
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To consider some possible decay channels for them, such as ttN, r]N, KH, i] b N and so on, 
we estimate these decays through heavy beauty meson exchanges by means of box diagrams 
as m 

Refs.flflfl- 

We neglect transitions to the hidden charm channels such as DT, C 
and -DA+, because they need t-channel exchange of too heavy vector meson constituted of 
charm and beauty quarks. The results for PB and corresponding VB channels are listed 
in Table |Vj Comparing results in Table [IV] and Table |Vj the influence of these additional 
coupled channels to the masses of predicted states is negligible. This is because the transition 
potential by exchanging heavy beauty vector meson is much smaller than the potential by 
exchanging light vector meson. 

TABLE IV: Pole positions zr and coupling constants g a for the states in (I, S) = (1/2, 0) sector. 



z R (MeV) 


9a 






BZ b 


BA b 


11052 


2.05 







B*£ 6 


B*A b 


11100 


2.02 






TABLE V: Mass (M), total width (T), and partial decay widths (T;) for (I, S) = (1/2, 0) sector. 



M (MeV) 


T (MeV) 




r, 


(MeV) 








ttN 


rjN 


rfN KT, 


Vb N 


11052 


1.38 


0.10 


0.21 


0.11 0.42 


0.52 






pN 


uN 


K*T, TN 




11100 


1.33 


0.09 


0.30 


0.39 0.51 





We also do not consider the coupled channel effect between VB and PB channels as in 



Ref.|l6|. The reason is that the transition potentials PB — > VB are much smaller than 
the potentials of PB — > PB or VB — > VB. Taking BH b — > B*H b through t-channel pion 
exchange as an example, the B*7rB coupling is proportional to (pb — Vi^^b* anc ^ ^ s zero 
in the static limit which ignores the three momenta of mesons and assumes e B * = (0, £*#.). 
Going beyond the static limit will give a non-zero transition potential but still much smaller 
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than its diagonal partners. This has been demonstrated by the production rate of J/ip/r) c 



from pp collisions in Ref. 



161 ] . The cross section for pp — > ppj/if) through DT, C bound state 



is smaller than that for pp — > ppi] c by more than an order of magnitude for similar excess 
energies. Therefore, the coupled channel effect between VB and PB channels is expected 
to have negligible influence on our predicted states. 

One problem associated to the beauty sector should be addressed here. As shown in 
FigfJl the loop functions of the hidden beauty sector, calculated with the cut-off or with 
dimensional regularization, show a quite different energy dependence and cannot be made 
similar over a reasonable range of values, as is the case for the hidden strange sector. This 
is due to that the on-shell momentum in the beauty channel shows a much stronger energy 
dependence than in the lightest channel. The results listed in Tables IIV|IVl are obtained in 
the dimensional regularization scheme, where the subtraction constant is adjusted to the 
value of the A = 0.8 GeV cut-off loop-function at threshold. However, the binding energy 
is found to be about 35 MeV for the BUb channel, which lies quite away of its threshold, 
where the real parts of the two loop functions are very different. This makes the choice 
of matching point for the two loop functions questionable. In order to get some feeling 
about the choice of the match point, it is also tried to match the two loop functions at 30 
MeV below threshold. Then the regularization subtraction constant moves from -3.715 to 
-3.774, and the binding energy moves from 35 MeV to 59 MeV. If the A = 0.8 GeV cut-off 
loop function is used directly, then the binding energy increases to 145 MeV. So the simple 
Valencia model for the beauty sector works not as good as for the hidden strange sector. 
The uncertainty for the concrete binding energies is quite large of the order of tens to a 
hundred MeV. But the qualitative conclusion for possible existence of bound state should 
be very solid. 

Then we discuss the (I, S) = (0, -1) sector. There are 3 channels, B s A b , BE b and BE' b . 
The masses of B, B s , E b and A b have been precisely measured and can be taken from Ref.jl|. 
v(Ib b = 5.366 GeV, itlb* = 5.4128 GeV and m= 6 = 5.7924 GeV. The has not been observed 



yet. Its mass has been predicted to be 5.922 GeV in Ref. [221] and 5.960 GeV in Ref. [23]. 
We choose a middle value 5.940 GeV in this paper. From Table [TTJ the BS b channel is 



decoupled from other two channels, so there should be a bound state for this channel, the 



same as corresponding vector-meson-baryon channel, For this channel, the results are 



listed in Table ED 
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For the coupled B s Ab and channels, the T matrix can be written as: 

rp _ 1 | V B s A b ->B3 b G BS b V BaAb ^. BSb 

with V = V BSb ^ BSb + Vl Ab ^ BSb G BsAb . 

The V is negative and hence provides an attractive potential. For = —3.71, one pole is 
found for the coupled-channel system, with mass between the two thresholds of B s At (10.986 
GeV) and BSb (11.071 GeV). The pole position depends on the value of as demonstrated 
in Table I VI I and can move to below the B s threshold when the magnitude of increases, 
such as for a M = —3.82 corresponding to the A = 1.1 GeV. 

TABLE VI: Pole positions z R with different for PB -)• PB in (I, S) = (0, -1) sector. 



A (GeV) 




z R (MeV) 






B s Ab and BEb 


BE' b 




0.7 


-3.68 


11030 - 0.60i 


11198 


0.8 


-3.71 


11021 - 0.59i 


11191 


0.9 


-3.75 


11004 - 0.49i 


11178 


1.0 


-3.78 


10990 - 0.24z 


11167 


1.1 


-3.82 


10970 


11151 



The coupling constants and the possible decay channels of these two resonances are listed 
in Tables IVTT1 and IVlIII for a M = —3.71. Similarly, the results for the corresponding vector- 
meson-baryon channels are also listed in Tables IVTll and IVIIII for a M = —3.71. 

Totally two A^* and four A* states are predicted to exist with masses above 11 GeV and 
very narrow widths of only a few MeV. The very narrow widths are due to the fact that all 
decays are tied to the necessity of the exchange of a heavy beauty vector meson because of 
hidden bb components involved in these states, and hence are suppressed. If these predicted 
narrow A^* and A* resonances with hidden beauty would be found, they definitely cannot 
be accommodated by quark models with three constituent quarks. Together with other 
possible A^* and A* states of other quantum numbers with hidden beauty, they should form 
a super-heavy island of the heaviest masses for excited nucleons A^* and excited hyperons 
A*. 
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TABLE VII: Pole positions zr and coupling constants g a for the states in (I, S) = (0, -1) sector 
for = —3.71. 

z R (MeV) g a 

B s A b BE b BE' b 

11021 -0.59a 0.14-0. Hi 2.27 + 0.004i 

11191 1.92 

B*A b B*E b B*E' b 

11069 -0.59i 0.14-0.12* 2.24 + 0.005i 

11238 1.89 



TABLE VIII: Mass (M), total width (T), and partial decay widths (Tj) for the states in (I, S) = 
(0, -1) sector for = —3.71. 

M (MeV) r (MeV) T; (MeV) 

KN ttE r?A if A KE rj b A B s A b 

11021 2.21 0.65 0.01 0.08 0.14 0.01 0.19 1.18 

11191 1.24 0.28 0.18 0.10 0.18 0.48 

K*N lo A cpA K*E TA B*A b 

11070 2.17 0.61 0.01 0.01 0.20 0.01 0.19 1.18 

11239 1.19 0.26 0.26 0.17 0.48 



IV. EFFECTS OF MOMENTUM DEPENDENT TERMS IN THE POTENTIAL 

For our model calculations in the last two sections, the static limit is assumed for the 
t-channel exchange of light vector mesons by neglecting momentum dependent terms as dis- 
cussed after the Eq.([T]). However, in Ref. 24j, dynamically generated open charmed baryons 
were studied by solving the Lippmann - Schwinger equation beyond the zero range approx- 
imation. The momentum dependent terms were found to have non-negligible effects on the 
results. In order to investigate the possible influence of the momentum dependent terms in 
this case, in this section, we use the conventional Schrodinger Equation approach to study 
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possible bound states for the channel by keeping the momentum dependent terms in 
the t-channel meson exchange potential. 

The deduction of the momentum dependent potential by the t-channel exchange of light 
vector mesons is straightforward. By keeping momentum dependent terms up to quadratic 
order with proper normalization factor for the Schro dinger Equation and including the 
vertex form factors as in the Bonn potential model 
is obtained as the following: 

r S C ah m\ 1 /Ay - 2 



251 ]. the effective S-wave potential 



V ab{Pl B^P 2 B 2 ) ~ 4/2 ^ + ^ [ A 2 v + ( p J 

m 2 p + 2m 2 B + AmprriB r 2 2m 2 B — m p 



, k + 7 2 7 ?)> ( 10 ) 

Ampm B V6mpm B J 

where k and q are defined as (p + p')/2 and p — p' with p and p' the initial and final momenta 
of the pseudo-scalar meson, respectively, in the center of mass system of the BHi channel. 
For simplicity, we assume the same cut-off parameter Ay for the p and u mesons. 

The effective potential for the Schrodinger Equation in the coordinate space, V(r), can 
be obtained by using the following Fourier-transformation formulae: 

A 2 -m 2 2 1 1 (e-' mT e~ Ar 2 2 e~ Ar 

A2 _ m 2 2 tS. 1 / -mr -Ar A<=~ Ar 

^{( 



A 2 + q 2 q 2 + m? Air \ r r 2 

A 2 -m 22 k 2 , 1 /m 2 e- mr A 2 e" Ar A 2 — m 2 ,Ar , e" Ar 



A 2 + g 2 g 2 + m 2 4i \ 4 r 4r 4 2 r 

1 , 2 e" mr e~ Ar A 2 -m 2 e" Ar 
'~ r 2 A~ 

Then we can solve the Schrodinger Equation 

fi 2 

(-— V 2 + F(r) - E)V{r) = 0, (11) 

— LI 

to find possible bound state with eigenvalue E and corresponding wave function \l/(f), and 
estimate the size of the system f with the formula 



r 2 drdVt^*(f)r 2 ^(r). (12) 

It is found that whether there exists a bound state depends on the cut-off parameter Ay. 
The results corresponding to various Ay values are listed in TablellXl 
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TABLE IX: Eigenvalue E and average size of system r vs the cut-off parameter Ay. 
Ay(MeV) 1100 1200 1300 1400 1500 1600 1700 1800 1900 2000 
E(MeV) - -0.85 -4.49 -10.5 -18.4 -27.9 -38.7 -50.5 -63.3 -78.9 
f(fm) - 2.36 1.19 0.86 0.70 0.60 0.53 0.48 0.44 0.41 

From the Table ITXj we can see that when the cut-off parameter Ay is 1200 MeV or larger, 
the effective potential can provide enough attraction to form a bound state. For Ay in the 
range of 1200 ~ 1800 MeV, the binding energy is in the range of 1 ~ 50 MeV with the 
average distance between two hadrons to be about 0.5 ~ 2 fm. The typical values for Ay 
in the Bonn potential are A p = 1400 MeV and A^ = 1500 MeV {25]. The binding energy 
corresponding to Ay = 1600 MeV is quite close to that obtained by the Valencia approach 
in the last section. This gives some justification of the simple Valencia approach although 
there could be an uncertainty of 10 ~ 20 MeV for the binding energy. 

According to Ref . . "the apparent radius of the pion as seen by the photon is de- 
termined almost completely by the intermediate p meson: the intrinsic pion size must be 
considerably smaller than the measured charge radius. In descriptions which explicitly in- 
clude the p meson, the pion can therefore be considered point-like for all practical purposes" . 
In our approach with the t-channel p meson exchange explicitly included, the D meson sim- 
ilar to the pion is expected to have very small size while the intrinsic radius of S c baryon is 



26]. With typical size r larger 



expected to be around 0.5 fm similar to that for the proton 
than 0.5 fm, our predicted hadron molecular state should not suffer much from internal 
structure of the constituents. 

For the BgAb-BIib coupled channel case, it is not so easy to use the Schrodinger Equation 
approach. Since the simple Valencia approach gives a consistent result for the BUb single 
channel case with Schrodinger Equation approach, we expect it also gives reasonable results 
for the B s A b -BE b coupled channel case. 

V. PRODUCTION OF N*~ b AND A* 5 IN pp AND ep COLLISIONS 

In order to look for these predicted super-heavy N*^ and A*^ states, we give an estimation 
of their production cross sections in the pp — > pprjb and ep — > epT reactions. The Feynman 
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diagrams are shown in Figj3j We also estimate the background of the pp — > ppi]b with N*^ 
replaced by the nucleon pole. 




FIG. 3: Feynman diagrams for the reaction pp — > ppgb and ep — > epT. 



The Lagrangians for the interaction vertices of these two reactions are as follows 27H29|: 



£nntt = gNNnN^T -ip^N + h.c, (13) 

£nn v „ = gNN^NjEip^N + h.c, (14) 

£n? % n-k = g N *_+ N7T N*rNr ■ + h.c, (15) 

£N;- b Nr lb = g^+^N^N^ + h.c, (16) 

£ ee7 = ierfej^eA^ + h.c, (17) 



em 2 „ 



A>7 = -j^ffA^ + h.c, (18) 
Jp 

C NkNp = gN^NpN^l^Ng^P^Jr ■ % + h.c, (19) 

C N *- b Nr = gN^NpN^lB^Ng^P^J^ + h.c. (20) 



pii pv 



with g^(P) = -g^ + 

In our model calculation, we only consider S-wave PB and VB interactions, so the spin- 
parity J p of our predicted N*^ for the PB channels is l/2~, and the NTt for the VB channels 
can be either l/2~ or 3/2~, but assumed to be l/2~ here for a simple estimation of rough 
production rate. The coupling constants of the Lagrangians can be either calculated from 
its corresponding partial decay widths or obtained from references. They are all listed 
in Table [X] For the NNrji, vertex, the width of rjf, has not been measured. Since both 
r] b and r] c couple to nucleon through two gluon exchange, we use the relation gNN^ ~ 
g NNric a A s {M r)b ) /a^(M Vc ) to estimate the gNNr) b with gNN^ determined from the decay width 
of i] c — > pp. 
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TABLE X: The coupling constants of involved vertices and corresponding widths used. 



Vertex 


T(MeV) 


Coupling Constant (g 2 /At:) 


ppir 




14.4 


bb F 


0.033 


1.03 x 10~ 5 


K^PVb 
bb r 


0.52 


1.81 x 10~ 3 


ee'y 




1/137 


IP 




2.7 [27] 


K-b + PP° 


0.030 


1.55 x 10~ 8 


KM 


0.51 


4.72 x 10~ 4 


PPVb 




1 x 10~ 6 



As usual, the off-shell form factors should be considered here. We use two kinds of form 
factors for mesons and baryons, respectively. 

(21) 



Fm — 


"A-M m M 


A 2 - r? ' 
1V M Pm 


Fn = 







4) 2 ' 



(22) 

m t 

where M stands for it or p, and A" stands for A"*^ or nucleon pole. Here Am = 1.3 GeV, 
A^ = 1.0 GeV. 

To produce the predicted A^(11052) in the pp collisions, the center-of-mass energy should 
be above 12 GeV. In FigJU the left figure shows our theoretical estimated total cross sec- 
tion for the pp — » pprjb reaction through the A"* 5 production vs the center-of-mass energy, 
with (dashed curve) and without (solid curve) including the off-shell form factors. As an 
estimation of background contribution to the N*^ production, we also calculate the corre- 
sponding cross section through the off-shell nucleon pole without including the form factors. 
The result is shown by the dotted curve. The contribution from the nucleon pole is much 
smaller than that from the N*^ production, because the nucleon pole is much more off-shell 
than N*^. The contribution of the nucleon pole with form factors becomes very small for 
the same reason, so it is not shown in FigJU This background reaction will not influence 
the observation of the N% production, especially for the energy range of 13 ~ 25 GeV. The 
cross section from A^ production is about 0.1 nb, which is much smaller than that for the 
corresponding reaction pp — > ppr\ c with N* 5 production [16| of about 0.1 fib. The main rea- 
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son is that both couplings of N*^Ntt and N^Nrjb are much smaller than the corresponding 
N* 5 Nn and N* 5 Nt] c couplings. These two vertices cause a reduction of about 2 orders of 
magnitude. In addition, because the center-of-mass energy here is much larger than that 
in the previous calculation for the r] c production, the propagator of exchanged rr further 
reduces the contribution. For the same reason, the contribution with form factors is much 
less than that without them. 



75 100 125 150 175 200 

,1 




100 125 150 175 200 

Mem. (GeV) 




20 22 24 

Mem. (GeV) 



FIG. 4: Total cross section vs invariant mass of system for pp —> ppr\^ reaction (left) and e~p — > 
e~pT reaction (right), with (dashed curves) and without (solid curves) including off-shell form 
factors, through production of the predicted N*^ resonances. The dotted curve is the background 
contribution from the nucleon pole for the pp — > ppr/f, reaction without including form factors. 



For the production of A^(lllOO) in ep collisions, the invariant mass of the system should 
be above 11 GeV. The right figure in FigJH shows our calculated total cross section for the 
e~p — > e~pT reaction vs the invariant mass of the system with (dashed curve) and without 
(solid curves) including form factors. The cross section of this reaction is much larger than 
that for the pp — > pprjt, reaction. The reason is due to the propagator of massless photon. 
The propagator of photon is given as the following: 

11 

pj ~ 2(m 2 e + piPjcosO - EiEf) ' (23) 

where the pt, Ei are the three-momentum and energy of initial e~, and pf, Ef for final e~. 9 
is the angle between initial and final e~ . When the directions of initial and final e~ are the 
same, i.e., cosO = 1, the value of Eq. fl23l) becomes very large because of the very small mass 
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of e~ . As the beam momentum of e~ becomes larger, the propagator of photon can reach 
very big value. For the invariant mass of the system less than 15 GeV, the cross section of 
e~p — > e~pT reaction is of the same order of magnitude as that of pp — > pprj b reaction. 



VI. SUMMARY 



In summary, the meson-baryon coupled channel unitary approach with the local hidden 
gauge formalism is extended to the hidden beauty sector. Two N*^ states and four A*^ states 
are predicted to be dynamically generated from coupled PB and VB channels with the same 
approach as for the hidden charm sector [3] . Because of the hidden bb components involved 
in these states, the masses of these states are all above 11 GeV while their widths are of only 
a few MeV, which should form part of the heaviest island for the quite stable N* and A* 
baryons. The nature of these states is similar as corresponding N* 5 and A* e states predicted 



in Ref. 



16] , which definitely cannot be accommodated by the conventional 3q quark models. 



Production cross sections of the predicted N*^ resonances in pp and ep collisions are 
estimated as a guide for the possible experimental search at relevant facilities in the future. 
For the pp — > pprjb reaction, the best center-of-mass energy for observing the predicted N*^ 
is 13 ~ 25 GeV, where the production cross section is about 0.01 nb. For the e~p — > e~pT 
reaction, when the center-of-mass energy is larger than 14 GeV, the production cross section 
should be larger than 0.1 nb. Nowadays, the luminosity for pp or ep collisions can reach 
10 33 cm _2 s _1 , this will produce more than 1000 events per day for the N*^ production. We 
expect future facilities, such as proposed electron- ion collider (EIC) [30j, to discover these 
very interesting super-heavy N* and A* with hidden beauty. 
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